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Griffiths Inequalities for Ising Spin Glasses on the Nishimori Line 
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The Griffiths inequalities for Ising spin glasses are proved on the Nishimori line with various 
bond randomness which includes Gaussian and ±J bond randomness. The proof for Ising 
systems with Gaussian bond randomness has already been carried out by Morita et al, which 
uses not only the gauge theory but also the properties of the Gaussian distribution, so that 
it cannot be directly applied to the systems with other bond randomness. The present proof 
essentially uses only the gauge theory, so that it does not depend on the detail properties 
of the probability distribution of random interactions. Thus, the results obtained from the 
inequalities for Ising systems with Gaussian bond randomness do also hold for those with 
various bond randomness, especially with ± J bond randomness. 

KEYWORDS: spin glass, Griffiths inequality, Ising model, Nishimori line, gauge theory, ther- 
modynamic limit 

1. Introduction 

The Griffiths inequalities make significant contributions to the understanding of phase 
transition for ferromagnetic Ising models.^) For one formulation, the Griffiths inequalities are 
written as 

dP 



and 

d{Sc 



> 0, (1) 

where Jb is a positive interaction and Sc = HiGC ^ii^i — ^'^) ^ product of Ising variables 
for arbitrary subset of the sites of the system. The above inequalities state that the pressure, 
P{= logZ), and the correlation function, {Sc), are monotonic increasing functions of the 
strength of any interaction, Jb- Using the two inequalities, various results can be proved, 
for example, the existence of the thermodynamic limit of the pressure per unit volume and 
correlation functions. They also give significant insights for the existence of ferromagnetic 
phase transition. 

For random Ising systems which have both ferromagnetic and antiferromagnetic interac- 
tions, inequalities analogous to the first Griffiths inequality have been proved with various 
bond randomness. ^"^-^ In these cases, Jb is a random variable, so that it has been proved 
that the pressure is a monotonic increasing function with respect to some parameter which 
controls the effect of an interaction term. Similar to the ferromagnetic cases, the inequalities 
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may, for example, be used for the proof of the existence of the thermodynamic hmit of the 
pressure density. 

As far as we know, however, the inequahty analogous to the second Griffiths inequality 
has only been proved on the Nishimori line, which is a restricted region in the phase diagram 
of random bond systems. Morita et al first proved the inequalities analogous to both the first 
and second Griffiths inequalities on the Nishimori line for Ising systems with Gaussian bond 
randomness. ^'^^ Using the two inequalities, the existence of the thermodynamic limit for the 
pressure per unit volume and correlation functions has been proved under various boundary 
conditions. Relations between the location of multicritical points for various lattices have also 
been derived. The proof, however, uses not only the gauge theory, but also the properties of 
the Gaussian distribution, so that the proof cannot be directly applied to the systems with 
other bond randomness, for example, it J bond randomness. 

Almost all the rigorous results obtained on the Nishimori line are proved, in essential, 
only using the properties of the gauge transformation, so that they hold not only for any 
lattice structure, dimension of the lattice, and range of the interaction, but also for various 
bond randomness. '''^^ For example, related to the second Griffiths inequality, the present 
author proved that the following inequality holds for finite itJ Ising systems, A and B, on 
the Nishimori line: 



proof was originally carried out for the it J Ising models, and used for deriving an inequality 
about the location of multicritical points of various lattices. However, the same argument can 
be applied to the systems with other bond randomness, since the proof only uses the gauge 
theory. 

Therefore, it is a natural question whether the existence of the two Griffiths inequalities on 
the Nishimori line is a special property of Ising systems with Gaussian bond randomness, or 
they hold without the detail property of the probability distribution of random interactions. 
The main results of the present paper is that the two Griffiths inequalities do hold on the 
Nishimori line for Ising systems with various bond randomness which includes both Gaussian 
and lb J bond randomness, where the proof essentially uses only the properties of the gauge 
transformation. Thus, the results obtained from the two inequalities for Ising systems with 
Gaussian bond randomness^'^-* do also hold for Ising systems with bond randomness considered 
in the present paper, especially, for the systems with it J bond randomness. 
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2. The Model and the Bond Randomness 

Let us first define several quantities. We treat the Ising spin system described by the 
Hamiltonian 

n = -Yl JaSa, (4) 

AcV 

where 

Sa = 1[ Si. (5) 

i<=A 

Here V is the set of sites, and the sum over A runs over all subsets of V among which 
interactions exist. The number of sites in A is arbitrary and may be different from subset to 
subset. The lattice structure is assumed to be reflected in the choice of A for which ^ 0. 
The partition function of this system is written as 

Z = ^exp(^/3AJA5A), (6) 

{sj Acn 

where we introduce local inverse temperature, Pa- The reason is that we use local inverse 
temperature. Pa, as a parameter which controls the effect of the interaction term, —JaSa, 
along the Nishimori line, as will be shown later. To investigate the property of a physical 
quantity at inverse temperature, /?, we must, of course, set all the local inverse temperature 
to Pa = P- The following argument, however, does not depend on the value of each local 
inverse temperature, Pa- 

Then the pressure of the system is written as 

P=[logZ], (7) 

where the configurational average over the distribution of bond randomness is written as [•••]. 

The probability distribution of the random interaction Ja is denoted by P{Ja)- In this 
paper, we investigate the case that the probability distribution, P{Ja), satisfies the following 
two conditions: 

P{-Ja) = P{Ja) eM-Wp,AJA) (8) 

and 

^^^ = {Ja-[Ja])P{Ja), (9) 

where Pp^A is a parameter which characterizes P{Ja)- 

Most of the probability distributions of random interactions investigated in the spin glass 
problems may satisfy the above two conditions. In the case of it J distribution with the ferro- 
magnetic bond concentration, pa{pa > 1/2), defining Pp^A as 

exp(2/3p,AJ) = 7^^, (J > 0) (10) 
J- ~ PA 
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we may rewrite P{Ja) as 

p{Ja) = paKJa-J) + {i-paWa + J) 

_ exp(/3p,^J) exp(-/3p,AJ) 



- Ur f{WA-J) + 5{JA + J))eMfiv,AJA)- 

For the Gaussian distribution with average, Jo,A) and variance, a a, defining Pp^A as 

^,A = ^ (12) 



^A 



we have 



1 f-i JA-JQ,A? 

1 (-^If^lA f-Jl 



1 z^l^ (-Jl\ J, 

71^ exp I 1 exp J exp(/3,^ J^). 



In both cases, we can easily see that P{Ja) satisfies two conditions, eqs. (8) and (9), by direct 
calculations. It is noted that, in the case of ±J distribution, [Ja\ = Jtanh(/3p^yiJ). 

In the present notation, the interaction term, —JaSa, satisfies the Nishimori condition 
when 

PpA = Pa. (14) 

Namely, when eq. (14) is satisfied for all the subsets {j4}, we may use the properties on 
the Nishimori line which have originally been proved by Nishimori, using the local gauge 
transformation/'^) 

Here, we show two properties of the probability distribution, P{Ja), which are often 
used in the rest of the paper. When the function, /(Ja), is an odd function of Ja, namely, 
/(-Ja) = -/(Ja), we obtain 

[/{Ja)] = [tanh(/3p,A Ja)/(Ja)]. (15) 

Also, when the function, /(Ja), is an even function of Ja, and monotonic increasing function 
of I I, we have 

[Ja tanh(/?A Ja)/( Ja)] > [Ja tanh(/3A Ja)] [/( Ja)] (16) 
The proofs of eq. (15) and ineq. (16) are shown in Appendices A and B, respectively. 

3. The inequalities on the Nishimori line 

Let us clarify the situation. We assume that all the interaction terms satisfy the Nishi- 
mori condition, eq. (14). In the following, we investigate the change of the pressure and the 
correlation functions with respect to arbitrary (Sb, and prove the following two inequalities: 

^[F1>0. (17) 
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and 

;^[{Sc)]>0, (18) 

where we denote the thermal average by angular brackets as (•••)• The above inequalities state 
that the configurational average of the pressure and the correlation functions are monotonic 
increasing functions of Pb on the Nishimori line. Here, we explain the role of Pb- When 
Pb is zero, there is no interaction term, —JbSb, in the system. Increasing the value of Pb 
makes the effect of the interaction term, —JbSb, larger, though we consider the restricted 
case that Pb always satisfies the Nishimori condition, Pb = Pp,B- Namely, increasing Pb for 
the present system corresponds to increasing the strength of an interaction for ferromagnetic 
Ising models. Thus, on the Nishimori line, ineqs. (17) and (18), play the same role as the two 
Griffiths inequalities do for ferromagnetic Ising models. 

Next, we briefly explain the procedures of the proof, since the proof of two inequalities 
can be carried out similarly. 

3. 1 The basic transformation of the physical quantity 
First, using the identity, 

exp{PBJBSB) = cosh(/3B'/B)(l + tanh(/3B ^5)5*5), (19) 

we rewrite the physical quantity, Q, so that the dependence of Pb may be seen explicitly. 
Corresponding to this procedure, we denote the configurational average over the distribution 
of bond randomness except Jb as [•••]'. The thermal average without the Boltzmann factor, 
bxp^PbJbSb), is also denoted as (• • • )'. 

3.2 The physical quantity on the Nishimori line 

Next, we rewrite the physical quantity, [Q], on the Nishimori line, using the local gauge 
transformation. Since we investigate the change of the physical quantity, [Q], in the situation 
where all the interaction terms including Jb always satisfy the Nishimori condition, we may 
use the properties on the Nishimori line, for example, [{Sc)] = [{Sc)'^], in any step of the 
proof. Actually, we perform the local gauge transformation for all the spin variables, and 
interactions except Jb- This can be done, since, in this step, there is no Boltzmann factor, 
^^^{PbJbSb)^ and the dependence with respect to Jb is explicitly seen, as will be shown 
later. 

3. 3 The change of the physical quantity on the Nishimori line 

Finally, we consider the change of the physical quantity, [Q], with respect to Pb- Using 
eq. (9), we rewite the total derivative of [Q] by Pb as 
d d 



J l-OO 

dR = IT / <^JbP{Jb)[Q\ 

dPB dPB J -co 

dJB [dp^P^^B)] [Q]' + J dJBP{JB) 



dp 



B 
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We evaluate the second term of rhs by direct calculations. For the first term of rhs, when Q 
is an even function of Jb and monotonic increasing function of | Js |, we have 

[{Jb - [Jb])Q] = [{Jb tanh(/?B Jb) - [Jb tanhi(3 b J b)])Q] > 0, (21) 

where we use eq. (15) and ineq. (16). It must be noted that, for the itJ Ising models, Q 
mentioned above becomes constant with respect to Jb, since = J^, so that eq. (21) takes 
the value, zero. 

4. The proof of the first inequality 

In this section, we prove the first inequality (17) following the procedures shown in the 
previous section. 

4-1 The basic transformation of the pressure 

Using the identity (19), we rewrite the pressure [P] as 



[P] = [logZ] = [logZ + [log(cosh(/3BJB))] + [log(l + tanh(/3BJB)(5B)') 

Here, Z' denotes the partition function without the Boltzmann factor, exp(/3B Jb-Sb). 
4-2 The pressure on the Nishimori line 

We expand the third term of rhs of eq.(22) as: 

n L 



(22) 



log{l + tanhipB Jb){Sb)')\ = ^^(-l)""'- {/3b J b){Sb 



(23) 



n=l 



On the Nishimori line, using the local gauge transformation and eq. (15), we can rewrite each 
odd term of the series as 



tanh2"-i(/3BJB)(5B) 



'{2n~l) 



tanh2"(/5BJB)(5B) 



'2n 



(24) 



where, the local gauge transformation is performed for all the spin variables and all the 
interactions except Jb- Then, we can further rewrite the third term of rhs of eq.(22) as 

1 



log{l +tanh{PB Jb){Sb) ) 



^ 2n(2n - 1) 

,71=1 ^ ' 

[HJb)], 



tanh^"(/3BJB)(SB: 



'2n 



(25) 



were we define F{Jb) as 
F{Jb) 



= + tanh(/3B JB)(5ij) ) log ( 1 + tanh(/3B Jb)(Sb) 

+ ^ ((l - tanh(/3B Jb)(5b)' ) log ( 1 - tanh(/3B Jb)(Sb) 



(26) 



By direct calculation, we have 



F{-Jb)=F{Jb). 



(27) 
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We also obtain 
d 



dJ. 



F{Jb) 



(3b{Se 



■log 



l + tanh{(3BJB){SB) 



> 0, {Jb > 0) (28) 



B 2 cosli'iPB Jb) \ l -tanhiPB Jb){Sb)\ 

regardless of the sign of {Sb)' ■ Namely, F{Jb) is an even function of Jb and monotonic 
increasing function of | | . 

^.3 The change of the pressure on the Nishimori line 

Now, we calculate the total derivative of the pressure, [P], by f3B on the Nishimori line. 
First, we have 



d 



^ -[log(cosh(/3BJB))] + 



df3B d^B d(3B 
The first term of rhs of eq. (29) is directly calculated as 
d 



log(l+tanh(/3BJB)(5B) 



(29) 



B 



-[log(cOsh(/3BJB))] 



= [{Jb - [Jb]) log(cosh(/3B Jb))] + [Jb tanh(/?B J^)] 
> 0, (30) 



where we use ineq. (21), since log(cosh(/?B J^)) is an even function of Jb, and monotonic 
increasing function of | |. 

Similarly, the second term of rhs of eq. (29) is calculated as 



d 



[{Jb-[Jb])F{Jb)] + 



log(l + tanh(/3BJB)(5B) 
Using, ineq. (21), it yields that 

[{Jb - [Jb)])F{Jb)] > 
For the second term of rhs of eq. (31), a direct calculation gives 



d 



B 



■F{J 



B 



d 



dl3_ 



B 



■F{J 



Bj 



Jb{Sb) 



log 



'l + tanh{PBJB){SBf 



>0, 



(31) 



(32) 



(33) 



2cosh2(/3BJB) \^ 1 - tanh{PBJB){SBy 
since it is easily seen that each term in the square brackets of rhs of eq. (33) is nonnegative 
regardless of the sign of Jb{Sb) ■ Thus, we obtain the first inequality (17). 

5. The proof of the second inequaUty 

We can similarly prove the second inequality (18), which states that the correlation func- 
tion is a monotonic increasing function ol (3b- 
5.1 The basic transformation of the correlation function 

Using the identity (19), we rewrite the correlation function as 

{Sc)' +iaMPBJB){SBSc)' 



^^^^ l + tanh(/?BJB)(5B)' 

In this subsection, for simplicity, we define xb and 

Xb = tanh(/3BJB), 



(34) 
(35) 
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We further rewrite the correlation function as 

(Sc)' + xbHSbSc)' - {Sb)'{Sc)') - xI{SbSc)'{Sb) 



{Sc) 



I-xHSbY' 
xbHSbSc)' - {SB)'{Sc)')-xl{{SBScy - 



(Sc) + 



I-xHSbP 

Then, the configurational average of the correlation function can be written as 

\{SbSc)' - {Sb)'{Sc)'){xb-xI{Sb)') 



■ (37) 



[{Sc)] = [{Sc) ] + 



[{Sc) ] + 



1-xI{Sb)'' 
{{SbSc)' - {Sb)' {Sc)'){xp^BXB - x%{Sb)') 



(38) 



where we use eq. (15) for the second term of rhs. 
5.2 The correlation function on the Nishimori line 

On the Nishimori line, since [ip^B = Pb, and Xp^B = xb, the correlation function can be 
rewritten as 

^xIHSbSc)' -{Sb)'{Sc)'){1-{Sb)') 



[{Sc)] = [{Sc) ] + 



(39) 



Performing the local gauge transformation for all the spin variables and interactions except 
Jb, we obtain 



xI{Sb) {S 



c) 



I-xHSbY' 



x%{SbSc) 
1-xI{SbP 

x%{SbSc)' {Sb 



xI{SbSc)^ 
1-xI{Sb)'\ 

x%{SbSc)' {Sb)' {Sc)' 



l-xl{SBr 



1-xI{Sb)'^ 



and 



xI{Sb)'HSc)' 



Thus, the second term of rhs of eq. (39) becomes 
xI{{SbSc)' - {Sb)' {Sc)'){l- {Sb)') 



xI{{SbSc)' -{SB)'{Sc)'f 



Hence, we obtain 

[{Sc)] = [{Sc)']' + 



ts.nh\pBJB){{SBSc) - {Sb) {Sc) )' 



(40) 
(41) 

(42) 

(43) 
(44) 



1 -tanh\PB Jb){Sb)'^ 

where, for the second term of rhs of eq. (44), it is easily seen that the term in the square 
brackets is an even function of Jb and monotonic increasing function of | Js | . 
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5.3 The change of the correlation function on the Nishimori line 

The total derivative of the correlation function by Pb on the Nishimori line can be written 

as 



d 



dPi 



c) 



{Jb - [Jb]) 



tanh\f3BJB){{SBSc)' - {Sb)'{Sc)')' 



+ 



l-tanh2(/3BJB)(5ij)'2 
d ( is.nh\l3BJB){{SBSc)' - {Sb)'{Sc)'?^ 



(45) 



dPB\ l-tanh2(/3BJB)(5B)'2 

Using ineq. (21), the first term of rhs of eq. (45) becomes nonnegative. For the second term 
of rhs of eq. (45), a direct calculation gives 

d ( is.nh\PBJB){{SBSc)' - {Sb)'{Sc)'?^ 



df3B\ l-tanh2(/3BJB)(5B)'2 

2JBt^n\i{l3BJB){{SBSc)' - {Sb)'{Sc)'? 
cosh2(/3B Jb)(1 - tanh2(/3B Jb)(Ss)'2)2 

Thus, we obtain the second inequality (18). 



> 



(46) 
(47) 



6. Summary and discussion 

We have proved two inequalities, ineqs. (17) and (18), for Ising spin glasses on the Nishi- 
mori line with various bond randomness which includes Gaussian and it J bond randomness, 
where the probability distribution of random interactions must satisfy two conditions, eqs. (8) 
and (9). The two inequalities, which correspond to the Griffiths inequalities for ferromagnetic 
Ising models, state that, along the Nishimori line, the pressure and the correlation functions 
are monotonic increasing functions of any (3b which controls the effect of the interaction term, 
—JbSb- The present results are an generalization of those by Morita et al for Ising systems 
with Gaussian bond randomness,^) where the proof uses not only the gauge theory but also 
the properties of the Gaussian distribution, so that it cannot be directly applied to the systems 
with other bond randomness. The present proof essentially uses only the gauge theory, so that 
it holds without the detail property of the probability distribution of random interactions. 

Using the present proof, the results obtained from the two inequalities for Ising system 
with Gaussian bond randomness^'^-* can be also derived for the systems with various bond 
randomness. In the research of the Ising spin glass problems, however, most studies have been 
carried out for the systems with Gaussian or it J bond randomness. Thus, we may insist that 
the most important physical consequence of the present paper is that it is found that the 
results obtained for Ising systems with Gaussian bond randomness^'^^ do also hold for Ising 
systems with it J bond randomness. 

Let us briefly explain several physical consequences on the Nishimori line for regular lat- 
tices which can be derived from the two inequalities, though they have already been explained 
by Morita et al for Ising systems with Gaussian bond randomness.^) 
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From the first inequality (17), we can show that the pressure has the well known super- 
additivity,^''^) namely 

i 

where [P]v denotes the pressure of the set of sites, V , and 

V = Y,Vi, (49) 

i 

since [P]v is obtained from 'Y^i[P]vi by adding random bonds among Vi and Vj. Thus, we can 
show the existence of the thermodynamic limit of the pressure density on the Nishimori line 
under free boundary conditions, assuming invariance by translation with respect to random 
interactions and the stability boundedness.'^'^-* For itJ Ising spin glasses with short range 
interactions, it is easily seen that the pressure has a definite stable boundedness. 

From the second inequality (18), we can show the existence of the thermodynamic limit of 
the correlation functions on the Nishimori line under free and fixed boundary conditions. For 
free boundary conditions, when we consider two finite sets of sites, V and V' so that V' D V , 
we get 

[{Sc)]v' > [{Sc)]v, (50) 

since V' is obtained from V by adding random bonds. Thus, we may say that the correlation 
function is a monotonic increasing function with the system size. With the fact that the cor- 
relation function is bounded by unity, we can assert the existence of the thermodynamic limit 
of the correlation function under free boundary conditions. We can also prove the existence 
of the thermodynamic limit of the correlation functions under fixed boundary conditions by 
similar procedure. (See Ref. 5 for details.) 

From the second inequality, we can also obtain the relation between the location of the 
multicritical points. When the lattice Li is obtained from the lattice L2 by adding random 
bonds, we get the following inequality for the magnetization 

[{s^)]l, > m]L„ (51) 

from which, we have 

TcLi>Tc,L2, (52) 

where Tc^LiiTc,L2) is the temperature of the multicritical point of the lattice, ^1(^2). 

Finally, we mention about the relation between ineq. (3) and the second ineq. (18) of 
the present paper including the work by Morita et al.^^ Inequality (3) only states that, for 
finite systems, the correlation function of some positive I3b is larger than that of f3B = 0, and 
cannot give the information between the values of the correlation functions of two different 
positive Pb- Compared to the above fact, from eq. (44), we can explicitly see how the value of 
the correlation funcion increases as Pb increases. Thus, ineq. (3) has clearly less information 
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than ineq. (18), which comes from the derivative of eq. (44) by Pb- However, ineq. (50), for 
example, can be derived using only ineq. (3), though it was not explicitly mentioned in Ref. 9. 
In the proof of ineq. (3), the condition for the probability distribution of random interactions 
is only one condition, namely, eq. (8), which states that the system has the Nishimori line 
itself. Thus, following the argument executed in Ref. 5, the existence of the thermodynamic 
limit of the correlation function on the Nishimori line under free boundary conditions may be 
proved for all the systems which have the Nishimori line in the phase diagram. The relation 
between the location of multicritical points of various lattices can also be derived for the same 
systems. 
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Appendix A: Derivation of eq. (15) 

In this appendix we explain the derivation of eq. (15). Through Appendices A snd B, we 
denote the configurational average over the distribution of bond randomness except Ja as 



Using eq. (8) and the fact that /{Ja) is an odd function of Ja, we have 

/"OO 

[expi-pp^AJA)fiJA)] = / dJAPiJA)exp{-Pp,AJA)[fiJA)]' 



-oo 
oo 



oo 
oo 



which implies 



Thus, we also obtain 



dJAP{-JA) exp(/3p,AJA)[/(JA)] 

~oo 

dJAP{JA) exp(-/?p,AJA)[/(-JA)]' 

dJAPiJA) exp(-/?p,AJA)[/(JA)]' 
= -[exp(-/3p,A^A)/(^A)], 

[exp(-/?p,AJA)/(JA)] = 0. 
exp(-/3p,A-^A) 



f{JA) 



0, 



(Al) 



(A-2) 



(A-3) 



_exp(/3p^A^A) + exp(-/3p,A Ja) " 
since (l/(exp(/3p^A<^A) + exp(— /3p_a<^a))/(<^a)) is also an odd function of J a- Therefore, it 
yields that 

exp(-/Jp,A^A) 



[f{JA)] = [/(Ja)]-2 



exp(/3p,A-^A) + exp(-/?p,A^A) 
[tanh(/3p,AJA)/(JA)] 



f{JA 



(A-4) 
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Appendix B: Derivation of ineq. (16) 

When two functions, /{Ja) and Jyi tanh(/3yi Jyi), are both even functions of Ja and mono- 
tonic increasing functions of | |, we have 

[{Ja tanh(/3^ Ja) - Jo tanh(/3A Jo))(/( Ja) - /(Jo))] > 0, (B-1) 
for any constant value, Jo, from which, it yields that 

[JAtanh{f3AJA)f{JA)] > Jotanh(/3AJo)[/(JA)] + ([Ja tanh(/3A Ja)] - Jq tanh(/3A Jo)) [/(Jo)]'- 

(B-2) 

Here, we can choose the value, Jo, so that it satisfies 

[Ja tanh(/3A Ja)] = Jo tanh(/3A Jo). (B-3) 
Substituting eq. (B.3) into eq.(B.2), we obtain ineq. (16). 



12/[T3] 



J. Phys. Soc. Jpn. Full Paper 

References 

1) R. B. Griffiths: J. Math. Phys. 8 (1967) 478, 484. 

2) F. Gucrra and F. TonineUi: Commun. Math. Phys. 230 (2002) 71. 

3) P. Contucci and S. Graffi: J. Stat. Phys. 115 (2004) 581. 

4) P. Contuccci and J. Lebowitz: Ann. Henri Poincare 8 (2007) 1461. 

5) S. Morita, H. Nishimori and P. Contucci: J. Phys. A: Math. Gen. 37 (2004) L203. 

6) P. Contucci, S. Morita and H. Nishimori: J. Stat. Phys. 122 (2006) 303. 

7) H. Nishimori: Prog. Thcor. Phys. 66 (1981) 1169. 

8) H. Nishimori: Statistical Physics of Spin Glasses and Information Processing: An Introduction, 
Oxford University Press (Oxford, 2001). 

9) H. Kitatani: J. Phys. Soc. Jpn. 63 (1994) 2070. 



13/[T3] 



